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Abstract 

We determine the anomaly free discrete R symmetries, consistent 
with the MSSM, that commute with SU(5) and suppress the \i pa- 
rameter and nucleon decay. We show that the order M of such Zjjy 
symmetries has to divide 24 and identify 5 viable symmetries. The 
simplest possibility is a lif symmetry which commutes with SO(10). 
We present a string-derived model with this T>f symmetry and the ex- 
act MSSM spectrum below the GUT scale; in this model 2^ originates 
from the Lorentz symmetry of compactified dimensions. We extend 
the discussion to include the singlet extensions of the MSSM and find 
and Z§ are the only possible symmetries capable of solving the 
\x problem in the NMSSM. We also show that a singlet extension of 
the MSSM based on a Z^k symmetry can provide a simultaneous so- 
lution to the [i and strong CP problem with the axion coupling in the 
favoured window. 



1 Introduction 



Supersymmetric extensions of the standard model (SM), such as the mini- 
mal supersymmetric extension, the MSSM, promise to eliminate the hierar- 
chy problem. However they also introduce serious potential problems and to 
be viable they must evade the /i-problem and the problem associated with 
new baryon- and lepton-number violating terms. This suggests that there 
should be an additional underlying symmetry capable of controlling these 
terms. Dangerous dimension four operators can be forbidden by R- or mat- 
ter parity [IH3], which is an anomaly free Z 2 subgroup of the continuous 
baryon minus lepton symmetry U(l)#_£. Dimension five proton decay oper- 
ators can be forbidden by 'baryon triality' [3], which combines with matter 
parity to 'proton hexality' |5JE]- The latter is the unique anomaly free dis- 
crete non-i? symmetry forbidding the dangerous operators while allowing the 
usual Yukawa couplings, the /x-term and the effective neutrino mass operator. 
Anomaly freedom is believed to be a necessary property of discrete symme- 
tries as otherwise quantum gravity effects may render them inefficient [IjlTH S]. 

However, there are two unpleasant properties of these traditional discrete 
symmetries. First, they do not allow to address the /x problem. Second, they 
do not commute with the symmetries of the grand unified theories (GUTs) 
SU(5) or SO(10) [10]. In [TTJ a discrete R symmetry was identified which can 
address the /x problem and commutes with SO (10). This Z^ symmetry is 
anomaly free through cancellation by the Green-Schwarz (GS) mechanism. 
In a recent paper [12] we have shown that this Tjf is the unique possibility 
which commutes with SO(10), and we have pointed out that it also solves 
the problem associated with dimension five proton decay operators. Further- 
more it contains matter parity as a Z 2 subgroup that is left unbroken after 
supersymmetry breaking. 

In this paper we extend the discussion to consider the possible discrete 
symmetries of the MSSM which commute with SU(5). As we shall see, there 
are only five possibilities with the simplest one being the Z^. Our analysis 
applies to singlet extensions of the MSSM as well. 

The paper is organized as follows. In Section [2] we prove that there are 
only five (generation independent) discrete Z^ symmetries which (i) com- 
mute with SU(5), (ii) allow the usual Yukawa couplings and dimension five 
neutrino mass operator and (iii) address the /x and proton decay problems of 
the MSSM. Section[3]is dedicated to a more detailed discussion of the simplest 
such symmetry, Zf. We present a globally consistent string compactification 
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with the exact MSSM spectrum below the compactification scale. The model 
exhibits the lif symmetry, which originates from the Lorentz group of com- 
pactified dimensions. In Section H] we discuss discrete R symmetries in singlet 
extensions of the MSSM. In a theory with the usual NMSSM couplings the 
discrete R symmetries can, apart from suppressing the proton decay rate, 
provide us with a solution to the NMSSM hierarchy problem. In a differ- 
ent singlet extension, in which the singlet couples quadratically to the Higgs 
bilinear, we will identify a unique discrete R symmetry capable of solving 
the \x and strong CP problems simultaneously. Finally, Section [5] contains 
our conclusions. In two appendices we present a re-derivation of discrete 
anomalies in the path integral approach and collect anomaly coefficients for 
discrete R and non-i? symmetries. 

2 Discrete symmetries of the MSSM 

In this section we discuss discrete symmetries of the MSSM which commute 
with SU(5) and can solve the fi problem. As we shall see, the assumption 
that matter 1am charges commute with SU(5) allows us to restrict possible 
Zm symmetries of the MSSM, as well as singlet extensions, to only few 
possibilities. We start in Section 12.11 by showing that one cannot address 
the /i problem with non-i? symmetries. In Section 12.21 we then turn to the 
discussion of discrete R symmetries, for which we prove that the order M 
has to divide 24. Finally, in Section \2.2.'2\ we classify all possible charge 
assignments. 

2.1 Non— R discrete symmetries 

We start by discussing non-i? symmetries. We show that such discrete sym- 
metries that are consistent with SO(10) or SU(5) relations for matter, i.e. 
universal charges for quarks and leptons, cannot forbid the /x term (cf. the 
similar discussion in [T3]). 

Consider a Z« symmetry under which the three generations of Q, U 
and E carry discrete charge gf while L and D carry where g labels the 
generation index. Our conventions are given in Appendix [B] If the Zm 
charges obey the even stronger SO(10) relations (i.e. q\ = the following 
discussion applies as well. The anomaly coefficients A 3 := ^4su(3) c -su(3) c -z A/) 
A2 '■= A SU (2) l -su(2) l -Zm' ^1 := Aj(i)y-u(i)y-z M and A := A grav _ grav _ ZM 



2 



are (cf. Equation ( 1B.17|) in Appendix 
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.4, 



A 



5 X> • <tfo + #) , ( 2 - la ) 

1 9=1 

5 £(3 • & + #) + £ (fe. + fcfc) > ( 2 - lb ) 

9=1 

3 

A i = \ E ( 3 • 9?o + (*) + | ^ ' + ' (2 ' 1C) 

9=1 

3 

A = J2{W-qi + ^-qQ+2q Hu +2q Hd , (2.1d) 

9=1 

where the sum runs over the generation indices g and and qn d denote 
the Jjm charges of the up-type and down-type Higgs doublets, respectively 
Anomaly freedom requires 

(Ai<i< 3 mod 77) = 24(^0 mod 77) = p (2.2) 



with p 7^ in the case of GS anomaly cancellation (cf. Equation ( 1A.23I) in 
Appendix IA.3I) . Here we define 

f M forModd, 
' ' \ M/2 for M even . K > 

Condition ( 12.2ft implies 

A 2 -A 3 = mod?/ (2.4) 

and hence, also in the case of generation-dependent Zm charges, 

2 + teJ = mod »7 • ( 2 - 5 ) 

On the other hand, the condition that the p term is allowed is 

q Hu + q Hd = modM. (2.6) 

We therefore see that, if we demand SU(5) relations for matter charges, a 
non-i? Zm symmetry cannot be used to address the p problem, even if we 
allow for GS cancellation of anomalies. 
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2.2 Discrete R— symmetries 

Having seen that non-i? symmetries cannot be used to address the /i problem, 
we turn to discuss discrete R symmetries. In this subsection, we derive 
constraints on the order M of Zf 7 symmetries that can solve the \x problem 
and accommodate the structure of the MSSM. 



2.2.1 A constraint on the order M 

After adding the contribution of the gauginos and gravitino the anomaly 
coefficients are 



A* 



A 



R 



^EMo + ?§)-3, 

9=1 

1 3 1 



9=1 
3 



9=1 



(QH u + qH d ) ~ 11 



(2.7a) 
(2.7b) 
(2.7c) 



A R 



-21 + 8 + 3 + 1 + ^ [10 (gf - 1) + 5 (<?§ - 1)] 

9=1 

t" 2 (q Hu + q Hd ~ 2) , 



(2.7d) 



where qio, %, qu u and qu d denote the R charges of the matter and Higgs 
superfields, i.e. matter fermions and Higgsinos have charges q — 1. 

In the case p ^ 0, the GS mechanism requires the presence of an axion, 
such that A R is to be amended by the axino/dilatino contribution (55 = —1). 

Subtracting the coefficients from each other leads to the universality con- 
ditions 



A* - Af 
A* -A* 



mod r\ rx qn u + qu d = 4 mod 2rj , 
mod 77 



rv 



3 
5 



(qH u + q Hd )-6 



mod r) . 



(2.8a) 
(2.8b) 



Equation ( 12.8a[) is equivalent to 



2 (to. + to 



2 + T]> 



(2.9) 



with an integer i. Inserting this into Equation ( 12.8bj) yields 
3 

- [i T) - 4] = k T) 



with another integer k. Altogether we find 



(2.10) 



[3£-5k] = 12/?7 = I 




(2.11) 



In both cases 24/ M has to be integer, i.e. M has to divide 24. Thus the 
possible values of M are 3,4,6,8,12 and 240 In what follows, we consider 
all these possibilities. 

2.2.2 Classification 

Given the constraints on the order M , it is straightforward to classify all 
phenomenologically attractive charge assignments. Here we assume that the 
charge assignments are family blind. Though not absolutely necessary it does 
ensure that the symmetry does not prevent mixing between families in the 
fermion mass matrix. The classification was done by a scan over all possible 
values of M. In addition to forbidding the /i term we require that 

1. Mixed gauge-Zfj anomalies cancel, i.e. ^4f<j< 3 = p mod 77; 

2. Yukawa couplings 1010 H u and 10 5 Hd as well as the neutrino mass 
Weinberg operator 5 H u 5 H u are allowed; 

3. i?-parity violating couplings are forbidden. 

Under these constraints the allowed charge assignments are given in Table |2~T1 

For completeness we note that there are only two more charge assignments 
that are allowed demanding just the first two conditions. They are given in 
Table £U 

One may ask whether there are additional discrete symmetries, such as 
t}^ 1 x 1j(p \ which cannot be written as single Tiff symmetries but also 
fulfill the three criteria above. The only candidates for such symmetries are 
based on the two patterns shown in Table 12.21 We find that by amending 

1 We exclude the case M = 2 since there are no meaningful order 2 discrete R symmetries 



(cf. e.g. [14]). 
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Table 2.1: Phenomenologically attractive charge assignments. The charges (f^ 
and q s ^ d are Higgs charges shifted in such a way that the anomaly coefficients Af 1 
(1 < i < 3) are manifestly universal, p is the universal value of the anomaly 
coefficients; p ^ indicates GS cancellation of anomalies. 
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28 
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Table 2.2: Charge assignments which satisfy only the first two criteria. Both 
assignments have p = 0. 

these assignments by the usual matter parity one arrives at the TL^ symmetry 
of Table 12.11 Hence our classification also comprises the TL { q ] x Z { p ] case. 
Of course, in extensions of the MSSM, extra states can enjoy additional 
symmetries. 

2.2.3 Dimension five nucleon decay operators 

Note that the third condition is sufficient to eliminate baryon and lepton 
number violation due to dimension four terms in the Lagrange density. How- 
ever in the MSSM at dimension five there are problematic operators allowed 
that generate nucleon decay. To be consistent with the bounds on nucleon 
decay these must be suppressed by a mass scale more than eight orders above 
the Planck scale, a major problem. However in the case of the 1^ symme- 
tries these operators are automatically absent. To see this note that the 
requirement that up- and down-type Yukawa couplings be allowed implies 

3<?io + <7s + 1h u + q Hd = 4 modM. (2.12) 

Combining this with Equation (I2.8ap gives 

3q 10 + q E = mod M , (2.13) 
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showing that (for M^2) the troublesome dimension five operators 10 10 10 5 
are automatically forbidden whenever the Yukawa couplings are allowed. 

2.2.4 The gravitational anomaly constraint 

For all charge assignments, the MSSM contribution to the gravitational 
anomaly is 

^(MSSM) = 7 mod//. (2.14) 

All cases except for M = 6 have p ^ and hence require the presence of an 
axion a. Call the multiplet containing the axion S, 

S\ =o = s + ia; (2.15) 

later we will identify S with the dilaton. From the coupling to the gauge 
fields J d 2 8 S W a W a one infers that the axino/dilatino has R charge —1. 
Therefore, after adding the axino/dilatino contribution we obtain 

A^(MSSM + axino/dilatino) = 6 mod r) . (2-16) 

The condition for anomaly freedom is 

^(A* mod ?/) = A? mod r, (2.17) 

for 1 < i < 3. Now, since Af e Z and since the order M, and therefore r], 
divides 24, this condition is equivalent to 

Aq = mod T) . (2.18) 

From Equation (I2.16P we see that the cases M = 4 and 12 are anomaly free. 
The case M = 6 is anomaly free with an axion that is singlet under Z^. All 
the other cases require additional states in order to cancel the gravitational 
anomaly. 

However this does not necessarily require additional states in the low 
energy spectrum. This is because states contributing to the anomaly can 
acquire mass when the symmetry is spontaneously broken. Since the R 
symmetry is broken in the hidden sector when supersymmetry is broken 
these states can acquire a mass of order the supersymmetry breaking scale 
in the hidden sector which can be as large as 10 13 GeV. With this in mind 
we will not consider the gravitational anomaly any further. 
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2.2.5 Compatibility with SO(10) 

By looking at the symmetries in Table 12.11 we observe that only the Z^ 
symmetry is compatible with a complete unification of quarks and leptons 
(Qio — Qs)- We will now show that also the other cases can potentially be in 
accordance with SO(10). The Zfj can be a mixture of the \J(l)x subgroup 
of SO(10) and an additional Zf f symmetry which commutes with SO(10). 

The crucial point is to realize that SO(10) has rank five, and SU(5) and 
the SM gauge group have rank four. Hence, there is an extra XJ(l)x factor. 
We will denote the XJ(l)x charge by Qx- The branching rules are [15] 

SO(10) D SU(5) x U(l) x , (2.19a) 
10 -)• 5 2 + 5_ 2 , (2.19b) 
16 10_i + 5 3 + l_ 5 . (2.19c) 

Consider an SO(10) GUT with an additional Zjy- symmetry as given in Ta- 
ble 12.31 If the SU(5) singlets contained in the 16 and 16 representations, ipu 
and ip-ff, attain VEVs, which have U(l)x charge Qx = ±5, we arrive at the 
following breaking pattern 

SO(10) x -> SU(5) x Zfj . (2.20) 

Note that there appears to be an additional Z5 symmetry, which does however 
not constrain any couplings since it is the non-trivial center of SU(5) (cf. 

mm)- 

In summary, we can obtain our Zf^ symmetries from an SO(10) GUT. 
However, the scenarios presented here are only toy models. First of all, 
further Higgs fields are needed to break SU(5) down to the SM. In addition, 
to obtain doublet-triplet splitting and get rid of dimension five operators 
larger Higgs representations are needed. Also, anomaly matching (cf. [16]) 
forces us to introduce extra representations because the value of p does not 
equal the one given in Table 12.11 

These considerations also show that the SU(5) relations are mandatory. 
Since SU(5) and the standard model gauge group have the same rank, there 
is no U(l) with which our Zf f s could mix upon breaking SU(5). 

3 A simple ZjJ symmetry in the MSSM 

In Table I2TT1 we survey all symmetries and charge assignments which commute 
with SU(5). The simplest one, the Zf , commutes also with SO(10). In what 
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(a) MSSM field con- 
tent. 



(b) Higgs sector. 
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(c) Ti^j— charges for different values of M, the 
resulting T^j— charge as a linear combination of 
the U(l)x charge Qx and the Zjy- charge r. In 
the last column we list the value of the SO(10)- 
SO(10)-S^ anomaly p. 
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Table 2.3: Field content of an SO(10) x Zf f model which can produce the sym- 
metries in Table [2~fl 

follows we will discuss this case in more detail. 

3.1 Non— perturbative terms 

The gauge invariant superpotential of the MSSM contains 

W = n H u H d + Ki LiH u 

+ Y l JH d U Ej + Y l d j H d Qi Dj + Yj>H u Q t U 3 

+ A% L i L i E k + \f k LiQjD k + X^l U i Dj D k 

+ ^fj 1 H u Li H u Lj + K^ ki Qi Qj Qk Li + nf^ M UiUj D k Eg 

+ K fjl Qi Qj Qk H d + K-[f k Qi U j E k H d + nf^ LiH u H u H d . (3.1) 
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We see immediately that the coefficients /i, K i} \[f k , X\f k , X\f k , K\f w nf^ 

K fjh K iji an d K f are forbidden by Zf perturbatively while Y^ du and 
are allowed. In what follows we will show that at the non-perturbative 
level /i as well as and K^ ke will be induced while the R parity violating 
couplings Ki and A as well as the k^ 3 ~ 5 ^ remain zero. The reason is that the 
latter are forbidden by a Z 2 subgroup of 7jf which is equivalent to matter 
parity. This subgroup is unbroken by the supersymmetry breaking sector 
and thus remains a symmetry of the full theory. 

Let us spell out the argument in somewhat more detail. Call the 7*^ 
transformation £, 

( : matter superfield — > i • matter superfield , 

Higgs superfield —> Higgs superfield , 

9 ->■ i-9, 

W -)• -W . (3.2) 

Now look at the transformation under which matter superfields trans- 
form with a minus, Higgs superfields go into themselves and 9 — > —9. The 
transformation fermion — > — fermion and 9 — > —9 is a symmetry of any SUSY 
theory, therefore ( 2 is equivalent to matter parity, and, in particular, anomaly 
free with p = 0. One can use the path integral (cf. Appendix |A|) to show 
that correlators that vanish due to a non-anomalous symmetry with p = 
also vanish at the quantum level. Therefore, the matter parity subgroup 
contained in the %f will not be violated by quantum effects. 

On the other hand, correlators which are only forbidden by but not 
by Z2, i.e. which are invariant under £ 2 , can be non-trivial at the quan- 
tum level. A convenient way to parametrize effective couplings describing 
these effects involve the S field, which shifts under the Zjf symmetry as (cf. 
Equation flA~22|) ) 

S -)• S+ 1 - A GS . (3.3) 
The discrete shift of S is given by (cf. Equation f ]A.23j) ) 

A GS = -L(A G _ G _ zf mod 2) = (3.4) 
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with v 6 Z. This allows us to write down terms 



= exp I — 8ir 



2 l + 2n 



S) B + JlH u H d + K, 




Qi Qj Qk Lg 



1 + 2u 




(3.5) 



with some coefficients B , ~p and Klj'jJ and n G Z. Such superpotential terms 
are Zf covariant, i.e. the exponential transforms with a minus under Zf while 
the terms in the square brackets are invariant. Due to the fact that S enters 

g?T 2 l + 2n 1 

the gauge kinetic function, these terms are proportional to e 1+2v ^ . For 
n = v = they can be interpreted as originating from t'Hooft instanton 
effects. The 8ir 2 in the exponential can also be obtained directly in a stringy 
computation [TSj. The crucial property of the non-perturbative couplings 
( 13. 5 p is that they are naturally suppressed. 

o 2 l+2n e 

The critical question concerns now the interpretation of the e 1 + 2v 
terms. So far we have shown that such terms are Z^ covariant. In the 
MSSM as a 'stand-alone' theory, SU(3)c or SU(2)l instantons can generate 
such terms, but their magnitude turns out to be very small. Whether or 
not further terms, with given n and v, appear depends on the model. Let 
us now make the very common assumption that there is a hidden sector 
that gets strong at some intermediate scale A. Then the non-perturbative 
terms related to the strong dynamics may well be the source of supersymme- 
try breakdown [THIEO]. Given non-renormalizable interactions between the 
MSSM and the hidden sector, communicated by some messenger fields, A sets 
the magnitude of the MSSM soft terms, m so f t ~ A 3 /M^, with M* being the 
messenger scale. In such settings, holomorphic, i.e. superpotential, terms 
can also be induced by higher-dimensional operators. That is, the AW np 
terms can appear with magnitude m so ft ~ A 3 /M^, but in principle they may 
also be absent if there are no higher-dimensional operators connecting the 
MSSM sector with the hidden sector exhibiting strong dynamics. In other 
words, if the MSSM fields are singlets under the hidden sector gauge inter- 
actions, there is, a priori, no guarantee that the AW np terms appear with 
reasonable size. If the scale of MSSM soft terms is related to some hidden 
sector strong dynamics, we expect the holomorphic terms also to appear, un- 
less there are additional symmetries beyond Z^ that forbid such couplings. 
Assuming that the dominant non-perturbative scale is related to supersym- 
metry breakdown we expect that the AW np terms are of the order of the soft 
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supersymmetry breaking terms. We will mainly focus on gravity mediation, 
where M* = Mp and these terms are of the order of the gravitino mass m 3 / 2 
(in Planck units). Below in Section [3.51 we will present an explicit string 
theory example in which the non-perturbative \i term is directly connected 
to m 3/2 . 

At this point let us mention that for the case of discrete R symmetries we 
disagree with statements made recently in [13], where it was claimed that, in 
the context of gravity mediation, R symmetries will be broken at the Planck 
scale and be therefore ineffective. The claim relies on the observation that 
there are fields with Planck scale VEVs that break the R symmetry. The 
derivation of this result relies on the inequality |(>^)| < hf r \F\ (cf. Equa- 
tion (9) in [13]) where f r is the i?-axion decay constant. This was derived for 
the case of continuous R symmetries by taking the limit of an infinitesimal 
transformation [T3]. For the case of discrete R symmetries the inequality 
is no longer true and there is no requirement that R-non singlets acquire 
Planck scale VEVs. In this case the R symmetry can be broken at a much 
lower scale. This is the case in the supergravity examples discussed here. 
In them the breaking of the R symmetry occurs non-perturbatively at an 
intermediate scale in a hidden sector and it is the superpotential VEV (W) 
rather than a field VEV that is the order parameter for R symmetry break- 
ing. Since the superpotential only appears at the non-perturbative level it 
is small. Also all other R symmetry breaking terms are small. This applies 
also to other schemes such as the one discussed in [21], where a small (W) 
is a consequence of an approximate R symmetry. Here the R symmetry is 
broken perturbatively, but again the order parameter, i.e. the superpotential 
VEV, is very small. In conclusion, R symmetries are a useful tool also, or 
in particular, in gravity mediation, where the same parameter, the small su- 
perpotential VEV, both sets the scale of soft masses and cancels the vacuum 
energy. In what follows, we discuss how the connection between the AW np 
terms and m 3 / 2 arises in the scheme of Kahler stabilization. 

3.2 Dilaton stabilization and supersymmetry breaking 

At the present stage of the discussion, the S field has no potential and su- 
persymmetry is unbroken. An economical way to rectify this situation is 
to invoke the stringy scheme of Kahler stabilization [22H25] (see also Ap- 
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pendix[C])o In this case the term of the form e~ bS represents a hidden sector 
gaugino condensate [2U] , which sets the scale for supersymmetry breakdown. 
According to the above discussion, in the presence of our symmetry 

o 1 + In . . 

6 = 8^—. (3.6) 

Let us discuss what that means in the case of a hidden SU(iV c ) theory with 
Nf chiral superfields in the N c + iV c representations. Here the coefficient b 
is given by 

3 3-Stt 2 

b = 2(3 = 3N^W f ■ (3 - 7) 
Therefore 

3 . . i±2n r3gl 
3N C - N f 1 + 2v ' y ' 1 

In the scheme under consideration, supersymmetry is broken by a non- 
trivial VEV of Fs- This leads to gaugino and soft scalar masses, following 
the pattern of the so-called "dilaton dominated scenario" [26] . This scenario 
has a number of phenomenologically attractive features. In particular, due 
to flavour universality in the soft breaking sector, it avoids the SUSY FCNC 
problem. Also, most of the physical CP phases, e.g. arg(v4*M), vanish which 
ameliorates the SUSY CP problem. However in the dilaton dominated case 
the vacuum structure may favour an unacceptable colour breaking minimum 
[27]. Other phenomenological aspects have been discussed in [28] . 

Moreover, the (non-perturbative) superpotential acquires a non-trivial 
VEV as well, 

(W) ~ e~ b{s) ± 0. (3.9) 

All gauge invariant terms which have been forbidden because they have zero 
R charge can now be obtained by multiplying them with (W). (W) will 
hence be the order parameter for R symmetry breaking. Inserting this in 
Equation (I3.5P we find that there will be a fi term of the order of (W), 
i.e. of the order of the gravitino mass m 3 / 2 , as well as K^ ki ~ 10~ 15 /Mp. 
On the other hand, terms which have odd 7jf charge cannot be obtained 
by multiplying them by e _fc ^; these are precisely the R parity violating 
couplings Ki, \(°\ and A^ 2 -* in Equation (13. ip . showing again that matter 
parity will not be broken. 



2 Alternatively, other stabilization schemes, such as racetrack mechanisms, may be ap- 
plicable here. 
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3.3 Phenomenology 

The suppression of the term leads to a situation in which dimension 
five proton decay will be unobservably small. Therefore, proton decay will 
proceed through dimension six operators mediated by gauge boson exchange. 

In settings with discrete R symmetries one should worry about the cos- 
mological domain wall problem [25]. The domain walls form at the stage of 
R symmetry breaking, typically the scale of supersymmetry breaking. For 
the case of gravity mediation this is at an intermediate scale of O(10 12 ) GeV. 
Provided the Hubble scale during inflation is below this scale, domain walls 
have sufficient time to form and then they will be inflated away. The require- 
ment that no domain walls are created after inflation translates in an upper 
bound on the reheat temperature Tr, which, given the other bounds on Tr 
in supersymmetric cosmology, appears rather mild. 

A discrete R symmetry may also be useful for inflationary scenarios. For 
example, in [30], it is argued that a Zf symmetry, with inflaton field (f> 
carrying R charge 2, can be used to guarantee that the inflaton potential is 
flat near the origin and give enough inflation^ 

In summary, for the case of gravity mediated supersymmetry breaking, 
non-perturbative effects naturally generate a /i parameter of the order of 
the gravitino mass. The symmetry ensures that the proton decay rate is 
well below the experimental limit and an exact matter parity is left that 
guarantees SUSY particles can only be pair produced and the lightest SUSY 
particle is stable. Thus one is left with the usual MSSM phenomenology with 
negligibly small corrections from higher dimension terms. 

3.4 literature 

A version of the Z^ symmetry has been prosed by Kurosawa et al. [31] where 
the traditional version of anomaly constraints was imposed, i.e. the possibility 
of GS anomaly cancellation has not been taken into account. This lead to 
a setting in which extra light charged states were required to cancel the 
anomaly. 

The Zf symmetry with GS anomaly cancellation has also been discussed 
by Babu et al. [II] . There are several aspects in which our analysis differs 
from or goes beyond [TTj : 

3 Note, their definition of the order of the discrete symmetry differs from ours. What 
they call Zf we call Zf . 
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1. We discuss the uniqueness for the first time. 

2. We point out that the Z4 symmetry also suppresses dimension five 
operators. 

3. We present the first discussion of non-perturbative !Zf violating effects. 

4. Related to the previous point, Babu et al. only discuss generation of the 
/1 term by the Giudice-Masiero mechanism. However, there will also be 
a holomorphic non-perturbative (i.e. Kim-Nilles type) contribution. 

5. In [TT] the is argued to originate from an 'anomalous' U(l)#. We 
are not aware of a model in which such an 'anomalous' U(l)# appears 
in string models. However, we also cannot rule out this possibility. 

6. We present a detailed discussion of the mixed hypercharge coefficient 
At. 

7. Babu et al. do not discuss the gravitational anomalies. 
3.5 String theory realization 

In the above discussion we argued that, if some hidden sector strong dynamics 
was responsible for supersymmetry breakdown, also a fi term of the right 
size will be induced by this dynamics. In order to render our discussion more 
specific, we will now discuss an explicit, globally consistent string-derived 
model. Such models have the important property that they are complete, 
i.e. unlike bottom-up (or 'local') models they cannot be 'amended' by some 
extra states or sectors. This allows us to clarify whether or not a reasonable 
fi term will appear. 

Making extensive use of the methods to determine the remnant sym- 
metries described in [17], we were able to find examples realising the Ti± 
discussed in Section [31 based on the string model derived in [32] and similar 
models, with the exact MSSM spectrum, a large top Yukawa coupling, a 
non-trivial hidden sector etc. In what follows, we present an explicit exam- 
ple. 

Consider the MSSM candidate model of [32]. It is obtained by the com- 
pactification of the Eg x Eg heterotic string on a Z2 x 7L<i orbifold with an 
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representation 
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representation 


label 
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(3,2;l,l,l)i 


Q 
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(3, 1; 1, 1, 1)_| 


U 
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(3,1; 1,1,1) i 


D 




5 


(3,l;l,l,l)_i 


D 


7 


(l,2;l,l,l)_i 


L 




4 


(1,2; l,l,l)i 


L 


3 


(l,l;l,l,l)i 


E 




33 


(1,1; 1,1,1)0 


N 


5 


(1,1;3,1,1) 


X 
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(1,1;3,1,1) 


X 


6 


(1,1; 1,1,2)0 


Y 




6 


(1,1;1,2,1) 


z 



Table 3.1: Spectrum of the orbifold model from [32]. The representations w.r.t. 
Csm x [SU(3) x SU(2) x SU(2)]hid, their multiplicities (#) and labels are listed. 



additional freely acting Z 2 . At the orbifold point (where the VEVs of all 
fields are set to zero) the E 8 x E 8 gauge group gets broken to 

G SM x [SU(3) x SU(2) x SU(2)] hid (3.10) 

times eight U(l) factors. One of them, denoted by U(l) anom , appears anoma- 
lous, i.e. trQanom = 180 7^ using the normalization |i an0 m| 2 = 15. Hence, a 
one-loop Fayet-Iliopoulos D-term gets induced. Furthermore, the massless 
spectrum includes three generations of quarks and leptons and is summarized 
in Table 13.11 More details on the model can be found in |33] . 

Next, we choose a vacuum configuration in which the SM singlets 

{(pi} = X$, X3, X4, X 5 ,Yi,Y 2 , Zi, Z 2 , 

N u N 2 , N 4 , N 7 , N w , N 15 , N 16 , N 17 , 

N 20 , N 21 , N 25 , N 27 , N 28 , N 30 , N 32 , N 33 } (3.11) 

attain VEVs. These fields are charged with respect to the hidden sector 
gauge group, the U(l) factors and several discrete symmetries (i.e. this or- 
bifold compactification provides three 71^ symmetries reflecting the discrete 
rotational symmetry of the three 1j 2 orbifold planes and six Z2 factors coming 
from the space group selection rule, see Appendix B of [32]). Hence, the 4>i 
VEVs of Equation (13. lip break these (gauge and discrete) symmetries and 
it turns out that 

G S m x Zf x Z 2 (3.12) 
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Table 3.2: x 1j 2 charges of the fields with SM quantum numbers. 



remains unbroken, where is a mixture of an orbifold l^f and other symme- 
tries]^ The x 7L 2 charges of the SM charged fields are listed in Table 13.21 
From there we see that the factor gives a stringy realisation of the 
symmetry described in a bottom-up approach in Section [31 Furthermore, 
the (pi VEVs also provide mass terms for the exotics, which are massless at 
the orbifold point, and allow us to cancel the Fayet-Iliopoulos D-term. 

In detail, we find four Fayet-Iliopoulos monomials, i.e. monomials that are 
gauge invariant except for a total negative U(l) anom charge such that their 
VEVs can cancel the positive Fayet-Iliopoulos term in the D-term potential 
(cf. [34] )• The monomials read 

{N&YM, N* 8 Z X Z 2 , iV 3 4 3 ViV 2 , N^Z X Z 2 ) , 

with Qanom(FI monomial) = — 15. Further, we find monomials with zero or 
positive U(l) anom charge involving all fa fields from Equation (13. lip . Hence, 
this represents a D-flat configuration. 

Matter fields are identified as fields with charge 1 , see Table I3.2[ and 
are given by Q x , Q 2 , Q 3 , U\, U 2 , U 3 , E u E 2 , E 3 , D 3 , D 7 , D 8 , L 2 , L 6 and 
L 7 . An inspection of the discrete charges of the Higgs candidates, i.e. the 
remaining L and L fields, reveals that there is one massless Higgs pair at 
the perturbative level. Unfortunately, the additional Z2, which we cannot 
break, forbids some Yukawa couplings such that the charged lepton and d- 
type Yukawa couplings Y e and Yd have rank 2. 

4 Meanwhile a very similar string model exhibiting vacua without the extra has been 
found [I!]. 
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More explicitly, we have computed the couplings using the well-known 
string selection rules [SHIES], which have been extended to the case of non- 
local GUT breaking [32] • For the Higgs mass matrix we find 
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11 
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#io 
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(3.13) 



where for example 3 denotes a sum of known monomials in the fields of 
Equation (13. lip starting at degree three. Therefore a linear combination of 
Li and L 2 as well as a linear combination of L\, L4 and L5 remains massless. 
The mass matrix of the extra colour triplets 



{6,} = {D^D^D^D^Ds} 



(3.14a) 
(3.14b) 



is 
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extra triplets 
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#20 
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13 












#28 








#10 


3 
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#33 
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#10 


) 



(3.15) 



So we see that most exotics decouple at the linear level in the VEV fields 
<pi, one pair of exotic triplets gets masses at order nine in the fields and 
another one at order three. One may speculate that this leads to the presence 
of colour triplets somewhat below the GUT scale, which may account for the 
fact that, within the MSSM, the strong fine structure constant = gf/(47r) 
turns out to be about 3% smaller than a\ and a.2 at Mqut- An important 
feature of the 7jf symmetry is that integrating out the triplets does not give 
rise to dimension 5 proton decay operators, as each triplet that couples to 
quarks and leptons pairs up with a triplet that does not (cf. the similar 
discussion in [37]). In other words, the mass partner 5 of a triplet 5 that 
couples to Qi L k (and therefore has charge 0) cannot couple to Qi Qj 
(FigureH). 
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Figure 1: Absence of induced QQQL operators, 
term is forbidden by . 



Either one vertex or the mass 



The Yukawa couplings are 
/ 1 

Y u = U 4 4 I +L 2 I 4 4 I , (3.16a) 
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Y d = L x 22 +L 4 1 (3.16b) 



Y^ ■ (3.16c) 



As already mentioned, and Y" e have rank 2. 

A crucial property of the string embedding is that, unlike in the bottom 
up approach, we have obtained an understanding of the origin of the 7L\ 
symmetry. In stringy language discrete R symmetries originate from what 
is called "^/-momentum conservation" [35], which reflects discrete rotational 
symmetries of compact space-time dimensions. In our orbifold we have three 
T 2 /^2 orbifold planes which can be rotated against each other by 180°. Each 
rotational symmetry manifests itself as a discrete order four R symmetry 
in the effective field-theoretic description of the model. The remnant 7jf 
discussed above is a linear combination of such symmetries and other discrete 
symmetries, either coming from the space group selection rule or emerging 
from continuous U(l) symmetries through spontaneous breaking (see the 
general discussion in [T7]). Note also that such discrete R symmetries can 
already appear anomalous at the orbifold point [38]. That is, no mixing with 
the so-called anomalous U(l) is required to obtain a symmetry whose 
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anomalies are canceled by the GS mechanising 

Perhaps the most important property of the model is that there is a 
proportionality between the holomorphic mass term connecting L\ and L\ 
and the superpotential VEV. This relation can be derived in an SU(6) orbifold 
GUT limit of the model, where it emerges due to gauge invariance in extra 
dimensions [39] (cf. also the field-theoretic discussion in [ffilUT]). Let us 
comment that the same SU(6) gauge symmetry also enforces the tree-level 
equality between the gauge and top- Yukawa couplings [32]. We hence see 
that, at least in this model, the superpotential VEV is both a measure for 
the gravitino mass, as usual, and the fj, term. 

4 Singlet extensions 

In Section 12.21 we have shown that the requirement of universality for the 
mixed gauge anomalies constrains the order M of a potential Zjy- symmetry 
to be a divisor of 24. As we have seen, this analysis carries over in an 
obvious way to singlet extensions of the MSSM, since additional SM singlet 
fields cannot change the constraints coming from the mixed gauge anomalies. 
In such extensions the MSSM subsector still has to obey the criteria derived 
in Section [221 However, the extra (singlet) fields can be subject to additional 
symmetries. 

In what follows we concentrate on two simple singlet extensions in which 
one or two singlet fields, respectively, are added. The first part of the discus- 
sion, Section 14.14 is on the so-called NMSSM, in which the singlet couples 
to the Higgs bilinear and there are cubic self-interactions. Section 14.21 is 
dedicated to a singlet extension of the MSSM which is capable of addressing 
the strong CP problem. 

4.1 NMSSM 

In the NMSSM, there is one additional singlet N with superpotential 

^ = ^ms = s°m + * N H u H d + kN 3 . (4.1) 
Let us now consider what this implies for the order M of a Zjj^ symmetry. 

5 The generality of the %f symmetry connected with a TLi orbifold suggests that 7L^ 
invariant models may also be obtained in other orbifold constructions such as or 
7Li x 7L 2 orbifolds of T 6 . 
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4.1.1 Constraints from NMSSM couplings 

There are three different classes of symmetries for which the iV 3 -term of 
Equation (14. ip implies different charges for the singlet N, i.e. 



M 
M 



mod 3 

1 mod 3 



no N 3 term possible 



M 



M = 2 mod 3 q N 
with q N the Zf f charge of N. 



3 

2M + 2 



mod M , 
mod M 



(4.2) 
(4.3) 

(4.4) 



M = 1 mod 3. Let us first consider the case M = 1 mod 3. The term 
XNH u H d together with Equation (I2.8ap then implies 



M + 2 



mod M + (4 mod 2rj) = 2 mod M 



M + 8 



mod M . 



(4.5) 



This equation has only one non-trivial solution for integer M, namely M — 4. 
Note that in this case qN = 2 mod 4 and a linear term in N is also allowed 
in the superpotential. Strictly speaking this is not the NMSSM but it is 
viable if the linear term is very small. We will discuss later why this may be 
natural. 

Following the analysis of Section |2~21 and using Equation (14. 3p . the unique 
charge assignment compatible with the Weinberg operator is shown in Ta- 
ble ED 
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<?10 


OS 








4 


1 


1 








2 



Table 4.1: Charge assignments for the symmetry. 



This is exactly the Tjf symmetry which we discussed in Section [3j We 
have seen that the mixed gauge anomaly coefficients of this symmetry satisfy 
the Green-Schwarz condition. Of course the singlet does not change these 
coefficients, so the analysis still applies. 
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M = 2 mod 3. Let us now consider the case M = 2 mod 3. The term 
XNH u H d together with Equation (I2.8ap then implies 

2M + 8 



mod M . 



(4.6) 



The solutions to this equation are M = 2, 8. As we have noted earlier there 
are no meaningful M = 2 R symmetries. The M = 8 case however is very 
interesting since, in this case, = 6 mod 8 and the linear term in N is 
forbidden. Following the analysis of Section |2~21 and using Equation ( 14. 4p . the 
unique charge assignment compatible with the Weinberg operator is shown 
in Table E2\ 
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Qh u 




qN 
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Table 4.2: Charge assignments for the Zg 2 symmetry. 



As the singlet does not contribute to mixed gauge anomalies, we know 
already from Table [27T1 that the lif symmetry has A^ <i<3 (MSSM) = p = 1. 



4.1.2 The hierarchy problem 

Searching for possible Zjy- symmetries in the context of the NMSSM we found 
that there are only two potential candidates: a and a Zf symmetry. 
The lif symmetry is actually a subgroup of the Z^ symmetry, hence both 
symmetries are closely related. While the commutes with SO(10) the 
Zf only commutes with SU(5). In both cases all dimension four and five 
baryon and lepton number violating operators are forbidden (except for the 
Weinberg operator), consistent with what we found in Section |2T2"1 

A potential problem with NMSSM models arises because SUSY breaking 
breaks the R symmetry and in radiative order a linear term in iV is generated 
in the superpotential. If the coefficient of this linear term is larger than the 
square of the electroweak scale it will lead to a large VEV for the singlet iV 
and therefore to a destabilization of the SUSY solution to the gauge hierarchy 
problem. This has been studied in detail by Abel jl3] who showed that the 
only dangerous operators that induce divergent tadpoles arise either from 
even terms in the superpotential or odd terms in the Kahler potential. He 
also showed that an R symmetry can avoid such terms because of the different 
R charges of the super- and Kahler potential (cf. also [44]). From the charge 
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assignments of Table |4~TI and Table I4~2l for the singlet N and the Higgs fields 
it is easy to show that the super- and Kahler- potentials actually do have 
exactly this structure in both the and the case and so in both cases 
radiative corrections do not destabilise the SUSY solution to the hierarchy 
problem. 

The main difference between and is that the former allows a linear 
term even at tree level. Does this mean that it is necessary to have the full 

symmetry when building the NMSSM? In effective theories, such as those 
describing the massless degrees of freedom in string compactifications, the 
superpotential starts with cubic terms in the fields and the linear term only 
appears through the coupling of the singlet field to fields acquiring VEVs. 
If the only (non-moduli) fields, <fr, with VEVs above the electroweak scale 
are in the hidden sector the coupling will be suppressed by messenger field 
masses, M*, which may be as large as the Planck scale. Allowing for trilinear 
couplings to messenger fields as well as trilinear couplings between messenger 
and hidden sector fields and assuming no additional symmetries, the lead- 
ing term in the superpotential after integrating out the messenger fields is 
iV0 4 1 Ml with the messenger scale M*. Taking Planck scale messengers, the 
constraint that this should not disturb the hierarchy is that (0) < y/M^Mp 
which is satisfied if the dominant VEV comes from the SUSY breaking sector. 
In this case it is sufficient to impose just the symmetry when building 
the NMSSM. 

The role of the SM singlets ^ with i?-charge 2 (such as N for the 
case Zf) has recently been discussed in the context of singlet (moduli) sta- 
bilization [18]. There it was found that for a superpotential with generic 
coefficients the number of singlets with i?-charge 2 should not exceed the 
number of fields <j)q} with i?-charge since otherwise the F-term conditions 
would overconstrain the system. Moreover, the ipf fields pair up with an 
equal number of fields. That is, for generic superpotential coefficients 
one might not expect to find vacua with an unbroken Zf symmetry and a 
massless singlet with i?-charge 2. However, it is quite conceivable that there 
are symmetries between the F-terms. In such a situation the i/j^ —<fio mass 
matrix won't have full rank such that one is effectively left with one (or more) 
singlet (s) with i?-charge 2. It will be interesting to see if this situation can 
be realized in string models in which there are additional symmetries, such 
as -D4 [151H6], relating the superpotential coefficients. 
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4.1.3 Non— perturbative effects 

Non-perturbative effects may also be important in determining the low en- 
ergy phenomenology. From Equation (13.51) we see that the superpotential 
has a term of the form 

AW np = B e~ bS (4.7) 

with a constant b. This parametrizes the non-perturbative effects discussed 
above, and may be interpreted as a hidden sector gaugino condensate. It 
provides the order parameter for local supersymmetry and generates the 
gravitino mass 

(A^np) (AA) 

~mT ~ ~M ~ m3/2 • ( ) 

AW np has .R-charge 2 (cf. the discussion in Section [3]) and similar non- 
perturbative effects can contribute to further terms in the superpotential. 
The crucial property of the non-perturbative couplings is that they are nat- 
urally suppressed. To parametrize these effects we denote by a superfield 
Y a non-perturbative term of the form given in Equation (14.71) (scaled by 
the factor Mp 2 ) carrying .R-charge 2 and we construct the superpotential 
involving Y that is consistent with the relevant R symmetry. 
The lowest superpotential terms in Y have the form 



AW%r = Y + Y 2 N + Y N 2 + Y H u H d 



~ m 3/2 M 2 + m 2 3/2 N + m 3/2 N 2 + m 3/2 H u H d , (4.9) 
AW%r = Y + Y 2 (N + Y N 2 + Y H u H d ) 



, m 3/2 , r2 , m l/2 



m 3/2 Mp + m 2 3/2 N+^N 2 + ^H u H d . (4.10) 



All of these terms have magnitude determined by the gravitino mass scale. 
For gauge mediation this scale can be very small and these terms negligible. 
For gravity mediation however the gravitino mass scale is the scale of super- 
symmetry breaking in the visible sector and the unsuppressed terms cannot 
be neglected. In this case, the magnitude of the AW%r terms is such as to 
reproduce the superpotential of the S-MSSM [EJUS], where, apart from the 
usual NMSSM couplings also holomorphic mass terms for the singlets and the 
Higgs fields of the order 7713/2 are introduced. This extension of the SM has 
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been shown to significantly reduce the fine tuning needed to accommodate 
the LEP Higgs mass bound jUJHH]. Our analysis yields a justification for 
the small holomorphic terms, which have so far just been imposed by hand. 

Interestingly the form of the non-perturbative effects is very sensitive 
to the underlying symmetry. For the case of lif there are additional un- 
suppressed linear and quadratic terms in N as well as a non-perturbative 
contribution to the Higgsino mass. For the case of only the linear term in 
iV is unsuppressed. Because the magnitude of all these terms is determined 
by the gravitino mass they will not disturb the SUSY solution to the hierar- 
chy problem. However, for the case of gravity mediation, the terms cannot 
be neglected and may be expected to significantly change the NMSSM phe- 
nomenology. Given the different non-perturbative terms appearing in the lif 
and we may expect these to have different phenomenological implications. 

4.2 Discrete R symmetries and the strong CP problem 

The axion solution to the strong CP problem remains the most convincing 
to date. Since it is based on the existence of a Peccei-Quinn (PQ) symmetry 
that forbids the \x term it is of interest to ask whether a discrete R symmetry 
can play the role of the PQ symmetry. Let us start by briefly discussing the 
role of axions in our setup with symmetries with particular focus on a 
possible solution of the strong CP problem. One potential candidate for such 
an axion is the universal Green-Schwarz axion, a = ImS, cf. Equation (12.15jl . 
As discussed in Appendix EH in the case of Kahler stabilisation we are — to 
leading order — left with a massless GS axion. However, as was e.g. shown 
in [19], the corresponding axion decay constant is of order the Planck scale 
if we demand the usual value for the unified gauge coupling. This is well 
outside the cosmologically allowed range 10 10 GeV < f a < 3 x 10 11 GeV, 
assuming no fine tuning in the initial axion VEV [49] . 

Kim and Nilles [50] have proposed an interesting model that naturally 
gives an axion decay constant in the favoured range. They achieve this by 
requiring that the coupling of the MSSM singlet field that contains the axion 
to the Higgs supermultiplets be quadratic with a superpotential of the form 

W = ^-N 2 H u H d . (4.11) 
M P v ; 

If the inverse mass scale, Mp, associated with this operator (the mediator 
scale) is taken to be the Planck scale with a = 0(1) an electroweak scale \x 
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term of (9(m 3 / 2 ) is generated if the singlet VEV ((N) = f a ) is in the desired 
range. This corresponds to (N) = O^m^^M-p) for the case of gravity me- 
diated supersymmetry breaking with 7713/2 close to the electroweak breaking 
scale. Furthermore, the theory has a global (accidental) U(l) PQ symmetry 
under which H u H d and the singlet N transform non-trivially. Hence, the 
VEV of N breaks the PQ symmetry and for the case of gravity mediated 
supersymmetry breaking the scale of breaking is of O(10 10 — 10 11 GeV). So 
the associated axion coming from the singlet iV has the right properties to 
solve the strong CP-problem |50j . 

One has to ensure that higher-dimensional operators, which explicitly 
break the PQ symmetry, do not spoil the solutionis An additional compli- 
cation in the construction of a viable model is that one needs an |7V| 6 term 
in the scalar potential to get an intermediate scale VEV of the correct mag- 
nitude if the soft supersymmetry breaking mass of the N field is 0(7713/2). 
Since a superpotential term N A breaks the PQ symmetry too strongly to give 
a viable axion at least one additional singlet (called X) is needed to generate 
the |iV| 6 term in the F-term potential. Let us consider whether the Zjjj- 
symmetries can give such a structure. 

In order to construct a viable model with intermediate scale breaking we 
start with the superpotential 

01 A72 TT TJ , P VAT3 



N 2 H u H d + —XN S . (4.12) 

M P u d M P K J 

which requires 1^ charges q^ = — 1 and qx = 5, cf. Table 12.11 Including 
soft SUSY breaking mass terms, this gives the potential 



r 



2a NH u H d + ^XN 2 



Mp " " M P 



2 



+ 



^ 3 

Mp 



2 



m 2 N \N\ 2 + m 2 x \X\ 2 
(4.13) 



that has the required \N\ 6 stabilising term. Provided m 2 x is negative, the 
field X acquires a VEV 

(X) ~ 0((N) 3 /M 2 ) , (4.14) 

where we have used (H u ) ~ (H d ) ~ ((N) 2 /M P ) = 0(m 3 / 2 ). The superpoten- 
tial in Equation f)4.12p has a global U(l) PQ symmetry under which H u H d , N 



6 A continuous R symmetry can protect the PQ symmetry at higher orders if it is broken 
only by the superpotential VEV and the intermediate scale VEVs of singlets [51] . 
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and X transform non-trivially. However the superpotential is not the most 
general one allowed by an underlying discrete R symmetry for this will allow 
additional terms of the form N p X q for p, q integer, where the values of p, q 
are constrained by the choice of Zjj^. Such terms will break the PQ symmetry 
generating a mass for the would-be axion. If the axion solution to the strong 
CP problem is to be maintained this contribution to the mass should be five 
orders of magnitude smaller than the corresponding contribution from QCD, 
Sm a < 10 -5 m^ CD < 10~ 9 eV [221121] • This puts a strong constraint on the 
discrete symmetry for it must be large enough to forbid the additional terms 
to a high order. 

Including the N p X q term in Equation (I4.12p one sees that it is the inter- 
ference between the last two terms in |-Fx| 2 that gives the dominant contri- 
bution because the VEV of X is smaller than the VEV of N. This term is of 
0{N p+3 X q ~ l /Ml +q+1 ) in the potential and gives a contribution to the axion 
mass given by 

5m 2 a = 0(N p+1 X q ~ 1 /MV +q - 2 ) = C(A^ p+3 ^ 2 /^ +3,? " 4 ) • ( 4 - 15 ) 

The leading terms for the candidate Zjy- symmetries (for which 5q — p — 2 = 
mod M) are N 2 , N 4 , X 2 , NX 3 and N 8 X 2 corresponding to p + 3q - 2 = 
0, 2, 4, 8 and 12 for M = 4, 6, 8, 12 and 24 respectively. The constraint 
5m 2 < 1CT 18 eV 2 is equivalent to p + 3q — 2 > 8 so we see that only is 
large enough to accommodate this method of simultaneously generating the 
/i term and solving the strong CP problem. 

Note that the singlets also induce baryon and lepton number violation 
as well as a small amount of R parity violation. In the case of the Z^ the 
leading order terms (cf. Table 12. ip are given by 



2 



XN 3 - , XN 

101010 5 and — r 1055. (4.16) 



With intermediate scale VEVs for N and (X) = 0((N) 3 /M$) the contribu- 
tion of these operators to nucleon decay is strongly suppressed lying below 
the irreducible dimension 6 operator contribution. Also the R parity viola- 
tion is negligible. In summary, a singlet extension of the MSSM with a 
symmetry provides us with a simultaneous solution to the \x and strong CP 
problems. The phenomenological implications for the Higgs structure may 
differ significantly from the NMSSM and remain to be analysed. 
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5 Conclusions 



We have discussed possible discrete symmetries for the MSSM which com- 
mute with SU(5). We have seen that, in order to address the fi problem, 
these have to be R symmetries. We have surveyed all possible discrete Zf f 
symmetries. Anomaly cancellation requires that the order M be a divisor of 
24. We identified 5 phenomenologically viable symmetries for the MSSM. 

The simplest of the 5 MSSM symmetries is a Zf which commutes with 
SO(10). This symmetry forbids all i?-parity violating couplings, dimen- 
sion five proton decay operators and the \x term at tree-level while allowing 
the usual Yukawa couplings and the neutrino mass operator. At the non- 
perturbative level the [i term and the dimension five proton decay operators 
are generated. We argued that in settings in which supersymmetry break- 
ing is related to some non-perturbative dynamics the \x term will be of the 
order of the MSSM soft terms. In particular, in gravity mediation we will 
have \x ~ fn 3 /2 and coefficients of the dimension five proton decay opera- 
tors ~ m 3/2/Mp, i.e. sufficiently suppressed. Thus the Zf symmetry 
provides us with a simultaneous solution to the arguably two most severe 
problems of the MSSM. We have discussed how to embed the Zf into string 
theory. Specifically, we have constructed a Z 2 x Z 2 orbifold with this Zf and 
the exact MSSM spectrum below the compactification scale, in which the 
lif originates from the Lorentz symmetry of compactified dimensions. At 
the non-perturbative level the /1 term is, due to the Z4 anomaly, generated. 
There is an exact matter parity and dimension five proton decay is well below 
experimental limits. 

We have discussed the role of discrete symmetries in singlet extensions of 
the MSSM. There are two possible symmetries consistent with the structure 
of the NMSSM, Zf and Zf, both of which are capable of solving the hierarchy 
problem. The Zf allows the usual couplings while forbidding the linear 
term for the singlet at the perturbative level. In the Zf case, one obtains 
holomorphic mass terms for the singlet and the Higgs at the non-perturbative 
level. We have argued that the size of such terms is of the order m 3 / 2 , leading 
to an S-MSSM-like scheme in which the smallness of the explicit mass terms 
for the singlets and Higgs finds an explanation. As another application we 
have discussed how discrete R symmetries can lead to approximate PQ U(l) 
symmetries capable of solving the strong CP problem. Given the upper 
bound on the order of Z^-, M < 24, we identified Z^ as the unique possibility 
for solving the \i- and strong CP problems simultaneously. 
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A Discrete anomalies in the path integral ap- 
proach 

In this appendix we re-derive Abelian discrete anomalies with the path in- 
tegral method, following [3S1I5I]- Among other things, we will describe how 
this allows us to understand the discrete version of the Green-Schwarz mech- 
anism. 

A.l Path integral derivation of anomalies 

Consider a theory described by a Lagrange density Jzf with a set of fermions 
\1/ = [ipW, . . . , ip( M '], where if)^ denotes a field transforming in the irre- 
ducible representation (irrep) of all internal symmetries. A general 
transformation ^ — > U ^ or, more explicitly, 



^(1) ■ 




(uw ^ 




■ ^(i) ' 




— » 












v UM J 




^(M) 



(A.l) 



which leaves =Sf invariant (up to a total derivative) denotes a classical sym- 
metry. 

A classical symmetry implies that certain correlators vanish at the clas- 
sical level. To see this, consider the correlator 

c n ,.. nM = ((^r---(^r<) . (a.2) 
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Now, if the field combination (ip^)™ 1 ■ ■ ■ (^( M )) nA / j s no ^ invariant under 
the symmetry transformation, we arrive at the (premature) conclusion that 

Cni...njvf 0. 

Classical chiral symmetries can be broken by quantum effects, i.e. have 
an anomaly. Specifically, consider a chiral transformation 

-> = exp (iaP L ) V(x) , (A.3) 

where a = a anom T anom with T anom denoting the generator of the transforma- 
tion and a anom being a parameter, and Pl is the left-chiral projector. 

We wish now to show that this implies vanishing correlators at the classi- 
cal level may appear at the quantum level. To this end, write the correlator 
as a path integral, 

C ni ... nM = Jv^VV(^r---(^ M Y M e iS , (A.4) 



where S denotes the action, which is left invariant under (1A.3|) . Now recall 



that under the transformation (1A.1I) the path integral measure undergoes a 
non-trivial change 



Vmvm -> J(a)VVVy , (A.5) 
where the Jacobian of the transformation is given by 

J (a) = expii [d 4 xA{a)\ . (A.6) 



The crucial observation is now that in the presence of a non-trivial Jaco- 
bian the full quantum correlator can be invariant. This is true regardless of 
whether the transformation I 1A.1I) is continuous or discrete, or whether it is 



gauged or global. 

The anomaly function A appearing in ( 1A.6|) decomposes into a gauge and 
a gravitational part [57H59] . 



A(a) = -4 g augc(a) + A grav {a) , (A. 7) 

with 

Aaugc(a) = — ^TrfaJ 7 .?] , (A.8) 

AUa) = -^^KKTr[a}. (A.9) 
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We have suppressed index contractions, i.e. = T^. Here Tp V = 
[Dp, D v ] is the field strength of the gauge symmetry, such that Tp V = (dpA u — 
d v Ap) for a U(l) symmetry, Tp V = F® v T a for non-Abelian gauge groups, 

and = 2£^ vpa J- p( j denotes its dual. Similarly, 1Z represents the Riemann 
curvature tensor and TZ71 = ^e^^lZ^'lZpax^ . The trace 'Tr' runs over all 
internal indices. 

Now we specialize to the case where a is a %m transformation. For the 
anomaly to be absent, i.e. J [a] = 1, we arrive at the conditions [381 151] 



G-G-Zm- Y,£(r (f) ) q (f) = mod 77, (A. 10a) 

r U) 

grav - grav - TL M : ^<7 (m) = mod r\ , (A. 10b) 

m 

where (cf. equation ( 12.31) ) 

M for M odd , 
M /2 for M even 

and q( m ^ denotes the Zm charge. The first sum runs over all irreducible rep- 
resentations r^> of G with Dynkin index £(r^) while the second sum runs 
over all fermions. Our conventions are such that £(N) = 1/2 for SU(Af) 
and £{N) = 1 for SO (AT). (lA.lOp are the traditional discrete anomaly con- 
ditions jU[S] with the difference that the T? M constraints do not appear; we 
will discuss 1? M anomalies separately in Section IA.4I 

A. 2 Green— Schwarz mechanism and re— derivation of 

#GS 

Consider a theory with simple gauge group G and an 'anomalous' Abelian 
gauge factor U(l) anom - Under U(l) anom the dilaton superfield S shifts ac- 
cording to 

S -> S+^6 GS A(x) (A.ll) 

with A denoting the U(l) anom transformation, i.e. the chiral superfields follow 
the rule 

$(/) g-iQinlmA $(/) _ (A. 12) 
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The corresponding transformation of the vector superfield V anom is 

Knom -> Knom + ^A-A^ (A. 13) 

with ReA|e =0 = a. In what follows, we derive the Green-Schwarz coefficient 
5qs from the requirement of invariance of the full action. 
The dilaton-dependent part of the Lagrange density is 



■^dilaton = — J ^d 4 # In (S + — 5 G S Kiiom) 



d 2 9 ^- Tr W a W a + h.c. 
4 

+ gravity terms . (A. 14) 

The first line of this Lagrange density is already invariant under the combined 
transformations (1A.11I) and ( 1A.13I) . The trace in the second line is supposed 
to run over all gauge factors, including U(l) anom . 

Decomposing the scalar component of the dilaton into a real and an 
imaginary (or axionic) part, 

S\e=o = s + ia , (A.15) 
leads to the usual couplings of the axion a 

& ^ - ^anom^anom " ^F'F* + ~11K , (A. 16) 

where F and -F anom denote the gauge field strength of G and U(l) anom re- 
spectively 

Hence, under a U(l) anom transformation with parameter a the axionic 
Lagrange density shifts by 

AJ^ axion = -^5 GS F anom F anom -^5 GS F a F a +^5 GS KK. (A. 17) 
lb lb 8 

The Green-Schwarz term Sqs can now be inferred by demanding that the 

transformation of the axion a cancels the anomalous variation of the path 

integral measure. The latter can be absorbed in a change of the Lagrange 

density 

A=^-anomaly F ansym F anom J 4.U(l)^ nom 



+ „ F a F a A G _ G _ V (i) 



a 

^^i^-^- ^grav-grav-U(l)anom ■ (A.18) 
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The coefficients A are the anomaly coefficients given by 

Aj(i)3 nom = ^E^™) 3 = ^ tr ^„o m , (A.19a) 

m 

^grav-grav-U(l)anom = ^ ] Qanam = ^ r Qanom i (A. 19b) 

m 

A G - G - Vi i Uom = E^ r(/) )^(L, (A.19c) 

r (/) 

where Qanom denotes the U(l) anom charge. The first two sums run over all 
left-handed Weyl fermions while the last sum runs over all irreducible rep- 
resentations r^' of G and £(r^) is the Dynkin index. 

The axion shift allows us to cancel the grav-grav-U(l) anom , U(l)^ nom and 
G - G - U(l) anom anomalies by demanding Aif anoma i y + A«^ axion = 0. This 
fixes the Green-Schwarz constant to be given by 



2n 2 5 GS = — trQ anom = -trQ^ om = A G _ G _ u(1)anom , (A.20) 



-trQ - 1 

2^ h Vanom g 

which is in agreement with the result obtained in a string computation [60] . 

A. 3 Discrete Green— Schwarz mechanism 

The Green-Schwarz mechanism also works if we replace U(l) anom by a dis- 
crete 7jm- In this case the parameter a is no longer continuous but a = 
with some integer n. Of course, there is no gauge field associated with the 
Zm, i.e. Equation ( 1A.13I) does not apply here. The discussion then goes as in 
the previous subsection. The discrete Green-Schwarz constant is now defined 
in such a way that under the 7Lm transformation of fields 

$(/) ^ e-^^/) (A.21) 
the dilaton shifts according to 

S S+'-Ags, (A.22) 

where Acs is fixed only modulo rj, 

7tMA gs = ^ -4 g rav-grav-z A/ = A G _ G ^ M mod 7} . (A.23) 
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The anomaly coefficients can be obtained from Equation (1A.19|) by replacing 
the U(l) anom charges Q^om by the Z M charges . Note that, unlike in the 
continuous case, the transformation of the axion is only fixed modulo 77. In 
the main body of the paper we obtain constraints on possible discrete sym- 
metries and charge assignments from the requirement that Equation (1A.23j) 
possesses a solution, i.e. that the Ag-g-tl m coefficients for different gauge 
factors G coincide modulo r\. 

A. 4 A comment on 7/? M anomalies 

If the discrete symmetry is embedded in a continuous symmetry, the univer- 
sality relations (1A.20j) also imply that 

A?rav-grav-Z M ~ 8 = mod X] . (A. 24) 

with 

= N ^(g (m) ) 3 , (A.25) 

m 

where M is a normalization factor compensating for the rescaling of the 
original U(l) charges to integer TLm discrete charges q( m \ However, this 
relation is an embedding constraint rather than a true anomaly constraint. 
That is, if this relation is not satisfied, this does not necessarily imply a 
non-trivial variation of the path integral measure [3H], arid therefore it does 
not mean that classically forbidden correlators will appear at the quantum 
level. 

The T? M anomaly constraints have lead to some confusion in the literature. 
Banks and Dine [9] gave an argument for why there are no T? u anomaly 
constraints. Following [S] they embedded the into a U(l) symmetry; 
however they broke this down to a lip.M symmetry and were able to show 
that, while there is only a Zm for the chiral states, the extra heavy states 
'see' a Ijp.m and can be chosen such that the anomaly conditions following 
from the U(l) constraints can be satisfied through their extra contributions. 

From this one might conclude that, in order to satisfy the 1? M anomaly, 
the true symmetry has to be TLp.u rather than %m- However, this is not nec- 
essarily the there can be discrete symmetries that cannot be obtained 
from continuous symmetries by spontaneous breaking in four dimensions. 
That is, the constraints from embedding Z« in non-anomalous continuous 
symmetries are sufficient to ensure anomaly freedom but not necessary. 
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In order to be specific, let us look at the stringy origin of the Tif sym- 
metry discussed in Section [3] of the main body of the paper. The 1^ has 
a clear geometric interpretation in terms of remnants of the Lorentz group 
of compactified dimensions. However, some of the states transforming non- 
trivially under our lif are twisted states. These states are chiral massless 
states which appear only after orbifolding, i.e. after we have broken the con- 
tinuous Lorentz symmetry down to a discrete subgroup. So it appears that 
there is no continuous interpolation between the continuous U(l) and the 
discrete Zf. Hence the derivation of discrete anomalies based on embedding 
discrete symmetries in continuous symmetries might not apply. On the other 
hand, the path integral method still works. This is consistent with the fact 
that in our orbifold construction there is no underlying 1i§ while we believe 
that the theory is UV complete. 



B 7Lm and anomaly coefficients 

We start by looking at the MSSM amended by ordinary, i.e. non-i?, discrete 
symmetries, where the fermions have the same charges as the superfields $^ 
and turn then to the discussion of discrete R symmetries. 



B.l Anomaly coefficients for non— R %m 

The anomaly coefficients for discrete non-i? Z m symmetries are well known [U 

Si, 

A G - G - %M = ^^')-? (/) . (B-la) 

r (f) 

A^-^-% u = ^q (m) ■ (B.lb) 
in 

These coefficients can be re-derived in the path integral approach [38] (cf. 
Appendix |A|). In Equation ( IB.laj) we sum over all irreducible representations 



of G while in Equation (IB. lb|) we sum over all fermions. £(r^) denotes 
the Dynkin index of the representation r^'. The discrete charges q are 
integers which are defined modulo M. Moreover, there might be mixed U(l) 
anomalies if the normalization of the U(l) factors is known. The coefficients 
are 

A u( i)-u(i)-^ = 5> (m) -(Q (m) ) 2 (B.lc) 
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with denoting the normalized U(l) charges. We will discuss this coeffi- 
cient in more detail below. 

Traditional anomaly freedom requires that for all anomaly coefficients 

A = mod 77. (B.2) 

However, discrete anomalies can be canceled by the Green-Schwarz mecha- 
nism, in which case one has to demand 

Ag-G-Z m = Aj(l)-U(l)-Z M = ^ Asrav-grav-ZM 

= p mod 77 (B.3) 

An important comment concerns the mixed U(1)-Zm anomaly coefficient 
( IB.lcj) . Mixed U(l) — U(l) — TLu anomalies are mostly ignored as they do 
not give meaningful constraints unless one knows the normalization of the 
charges [6j[6T]. Typically the sum in Equation ( IB.lcj) is not invariant under 



shifting some discrete charges by M. To see this, let us consider the example 
of hypercharge. We will denote the unnormalized U(l)y charge by Qy . The 
anomaly coefficient reads 

Ai = I (^y m) ) V m) = P mod 77 . (B.4) 

m 

We have the freedom to shift the Zm charges by integer multiples of M, i.e. 
we can define new Z M charges q' (m) = q {m) + fc (m) M with k {m) G Z. With 
the new charges the condition for anomaly freedom is 



Yl (Q^) 2 (q {m) + k {m) M) = p mod// (B.5) 

m 

A i + 5 M E fc(m) (^ m) ) 2 = P mod V- (B.6) 



We can choose the such that n is an arbitrary integer because, for 
example, Qy{E) = 1. Hence, we arrive at 

3 

A 1 = p nM + mr), meZ. (B.7) 

5 
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This can be rewritten as 

M odd : 5Ai = 5p + (5m - 3n) M , (B.8) 

M 

M even : 5A 1 = 5p + (5m - 6n) — . (B.9) 

Since 5m — 3n and 5m — 6n are arbitrary integers, we get 

5Ai = 5p mod rj . (B.10) 

B.2 Anomaly coefficients for 7L^ symmetries 

Now consider a Zjy- symmetry, under which, by convention, the superpoten- 
tial transforms as 

W e 27ri<? ^ /M >r (B.ll) 
with q^ = 2. Accordingly, the superspace coordinates transform as 

g _^ e 2»ri/Jlf 0j ( B12 ) 

such that d 2 # transforms oppositely to Superfields $^ = (j)^+y/2 9^ ^+ 
follow the law 

Correspondingly, the fermions transform as 

4>W = e 2 -(^-D/M^(/)_ (B14) 

For discrete R symmetries, the anomaly coefficients read (cf. Appendix IA~|) 

Aa-a-zfl = ^^(r (/) )-(g (/) -l)+^(adjG), (B.15a) 

r(/) 

V)-u(D-aa = E^) 2 -^'- 1 )' ( R15b ) 

m 

^ gr a V - gr av-^ = -21 + ^ dim(adj G) + #(U(1)) + ^(g( m ) - 1) . 

G m 

(B.15c) 

Here denote the 7L% charges of the superfields, the charges of the cor- 
responding fermions are shifted by one unit, q^f) = q^ — 1. In Equa- 
tion (IB. 15a|) £(adj G) = C 2 (G) represents the contribution from the gauginos, 
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#(U(1)) denotes the number of U(l) gauginos. The first and second term 
on the right-hand side of Equation flB.15cp represent the contributions from 
the gravitino and gauginos. A necessary condition for anomaly cancellation 
is the universality 



-4, 



G — G — S 



Aj(l)-U(l)-0« 

p mod T) . 



±a 

24 grav-grav-Z^ 



(B.16) 



p is a constant, which is related to the discrete shift (1A.22|) of the axion via 
p = irMA GS . 



B.3 Summary of anomaly coefficients 

The anomaly coefficients are given by 



.4 



IB) 



G-G-% 



Aj(l)-U(l)-^ 5 



A 

grav— grav— TLX, 



r (f) 



R) + £(adj G) ■ R 
-R), 



R 



-21 + ^dim(adjG) + #(U(l)) 

G 



(B.17a) 
(B.17b) 



(B.17c) 



where we distinguish between discrete non-i? (R = 0) and R (R = 1) symme- 
tries. #(U(1)) denotes the number of U(l) gauginos. As discussed above, the 
mixed U(l) — U(l) — anomaly is only meaningful if one knows the nor- 
malization. In general, the coefficient A ( R ) is not invariant under 

V(l)—U\L)—£ M 

( R) 

shifts of the 7L M charges by integer multiples of M. 



C A comment on Kahler stabilization 

A possible way to stabilize the dilaton is through non-perturbative correc- 
tions to the Kahler potential [22~||2"3] . Such corrections are expected to vanish 
in the limit of zero coupling and also to all orders in perturbation theory. 
The form of these corrections has been studied in the literature p24ll25| l62]. 



38 



For a favourable choice of the parameters, this correction may allow one to 
stabilize the dilaton at a realistic value, Re S ~ 2, while breaking supersym- 
metry [2H[251E2Hni]- A common parametrization of the non-perturbative 
corrections reads 

e K = e K 0+Q K np ^ ( C1 ) 

with K = — ln(2x), x = Re S, and parameters subject to K" > and 
p,q > 0. Supersymmetry is broken spontaneously by the F term of the 
dilaton, 

~Mp" ' (C - 3) 
For a single gaugino condensate, one has 

W = d exp f J , (C.4) 



2/3 

where 3/(2/3) = 8vr 2 /iV and d = -iV/(327r 2 e) for a condensing SU(iV) group 
with no matter. Note that the scalar potential is independent of Im S. That 
is, we are left with a GS axion. 

The problem with this scheme is that the vacuum energy at the local 
minimum is typically positive (cf. [63 1). Although Equation ( IC.2I) represents 
the 'standard' choice of the Kahler potential, there are no arguments that 
forbid additional terms of the same structure. That is, following Shenker's 
arguments (22] one may replace Equation flC.2j) by 

e K »* = ( Cl x pi/2 + c 2 x P2/2 ) e~ q ^ . (C.5) 

In fact there seems to be no reason for not writing even more terms in the 
parentheses. One can then tune the vacuum energy in the local minimum to 
zero by carefully adjusting the coefficients (Figure [2]). 

Of course there is still no reason for why the vacuum energy should vanish 
at the local minimum, but the above arguments may show that, in principle, 
the vacuum energy can be tuned to zero in this scheme. 
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Figure 2: Dilaton potential for {c 1 ,c 2 ,Pi,p 2 ,q} = {-28.8292,22.6129,2,3,4}. 
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